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Abstract 

We propose to use a 13 KeV antiproton beam passing through a dense cloud of 
positronium (Ps) atoms to produce an "beam". These ions can be slowed down 
and captured by a trap. The process involves two reactions with large cross sections 
under the same experimental conditions. These reactions are the interaction of p 
with Ps to produce H and the capture by H reacting on Ps to produce H . 
Once decelerated with an electrostatic field and captured in a trap the ions 
could be cooled and the e+ removed with a laser to perform a measurement of the 
gravitational acceleration of neutral antimatter in the gravity field of the Earth. 
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1 Introduction 



The measurement of the gravitational mass of antiparticles motivates physi- 
cists since over three decades. It is interesting to note that general relativ- 
ity would not contradict antiparticles to "fall up" in the gravity field of the 
earth [1]. Experiments to test such ideas have been proposed on positrons and 
antiprotons [2] but never succeeded. Several experiments at CERN are now 
producing neutral antimatter in the form of antihydrogen or antiprotonic he- 
lium atoms [3,4] for CPT tests. The trapping of neutral antihydrogen atoms 
is the next step for these experiments. Some proposals to measure the fall of 
these atoms have been presented [5]. The possibility to measure the free fall 
of positronium has also been studied [6] . 
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Fig. 1. Atoms/ions production scheme. 



Recently it has been proposed [7,8] to measure the gravity acceleration of an- 
timatter using the ion. This ion has the advantage that it can be cooled 
down to 13 /xK, a temperature suitable for a gravity experiment. The authors 
do not investigate in detail the production scheme but suggest the interac- 
tion of Ps Rydberg atoms with H atoms [9,10]. In this article we concentrate 
on the H production scheme and propose a "beam to cloud" experimen- 
tal configuration instead of a trapped system. This configuration allows the 
isolation of the from the production region. 

The positronium cloud hereafter referred to as "the target" is not contained: 
it is made of positronium atoms emitted from an Aluminium crystal surface 
which is bombarded by a fiux of positrons. The atoms are emitted in a direc- 
tion normal to the surface with a few mrd angular spread [11]. Their energy 
is 2 eV and their speed is ^1 mm/ns. The density of the Ps target is 
proportional to the positron fiux. As explained further, an amount of order 
10^^ positrons is needed in order to get the required density for the Ps tar- 
get. However, the short Ps lifetime requires this amount to be delivered in a 
few nanoseconds while the highest foreseen rates of slow positrons are at best 
10^^ s~^. The positrons have then to be accumulated in the crystal vicinity 
and accelerated toward the crystal. To counteract the effects of space charge, 
the positrons are accumulated into a small neutral e~e^ "plasma" for a short 
time. The plasma size is ~ 3mm^ x 1cm. Then an electrostatic field acceler- 
ates the positrons toward an Aluminium crystal where they are converted into 
Positronium atoms. 

A 13 keV antiproton beam is guided parallel to the crystal at a distance of 
^ 150 /im. The target length is crossed in 6 ns by the 13 keV antiprotons. 

The fiux of positronium atoms is maintained as long as the fiux of positrons 
toward the crystal is kept. It is the fiux of positronium atoms which constitutes 
the "target" . The overall layout for such a gravity experiment is presented in 
figure 1. 
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Aside from the gravity experiment, the technique to accumulate positrons 
in the plasma volume described in this article could prove useful as a first 
step toward the realization of a positronium [12] BEC, a 511 KeV laser [13] 
and the observation of the antimatter molecule Ps2. When Greaves-Surko 
traps [14,15,16] are available with a capability to hold 10^^ positrons and pro- 
vided that the neutral plasma can be held during the time required to empty 
the trap, one could produce a density of positronium of ^ 0.3 10^^ cm~^. An 
experiment to observe the stimulated annihilation process would then become 
feasible. 

If traps with such capability are not available, an alternative path toward a 
stimulated annihilation observation could be to use the neutral plasma in a 
MCEO trap [17]. The positronium atoms produced by the 3 body reaction 
are in Rydberg states. It takes a few /iS before they reach the ground state. 
There, the magnetic field couples the triplet state to the singlet state, which 
has a lifetime of 0.125 ns. Therefore the life time of the positronium produced 
inside the plasma is dominated by the decay time towards the ground state. 
Using an infrared laser onto the P* gas one can re-ionize some atoms while 
preventing the decay into states with a low level of excitation. Hence the life 
time can be increased and, in principle, longer trains of positrons could be 
filled into the plasma volume. 

This configuration is also well suited for spectroscopy experiments where 
atoms containing a positron, or a positronium atom [18,10] are produced by 
the interaction of atoms/ions with the target. The incoming antiproton beam 
is then replaced by an ion or atom beam. 

The path described in this article to produce an antimatter ion beam is the 
following: 

• the accumulation of positrons in a neutral e~e+ plasma, 

• the separation of e~ and by an electrostatic field and the interaction of 
the with an Aluminium crystal to produce a positronium cloud referred 
as the target, 

• the interaction of 13 keV (anti)protons with the target to produce the atoms 
and the ions as a beam. 

The beam to cloud configuration presented in this article is made of several 
devices which were developped separately for various applications: 

• the 10 MeV e~ beam on a thin foil to produce an intense source, 

• the buffer gas and Greaves-Surko trap, 

• the MCEO trap, 

• the Charge Focusing Aluminium Converter (CFAC). 

The present study takes the parameters of these devices as they appear in 
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the litterature and shows that with httle modifications, these devices can be 
assembled to produce an ion beam. 

The steps involved in the target production are discussed in the next section. 
The atom and ion production rates are presented in section 3. Several technical 
features are discussed in section 4. 



2 The positronium target 

The positronium target is obtained by accelerating positrons from an e~e"^ 
plasma toward an Aluminium crystal. The positrons hit the crystal with a 
kinetic energy above 40 eV to avoid elastic and specular reflection [19]. 

In the following we describe a way to produce the neutral plasma and the 
subsequent extraction and focalisation of the positrons from the plasma onto 
the Al crystal. 

We foresee two modes of operation for the creation of the neutral plasma: 

• the slow loading mode where the positrons are extracted continuously from 
the buffer gas section of the Greaves-Surko trap which cools them to room 
temperature 25 meV). The loading time is a few seconds. 

• the fast loading mode where the positrons are stored and cooled to 2 meV 
in a Greaves-Surko trap [16] and extracted in ^ 10 /iS. This fast extraction 
heats the positron beam. This mode of operation requires less confinement 
time for the neutral plasma and may allow higher plasma densities. 

The electrons are provided by a buffer gas or magnetic trap depending on the 
desired temperature. Since electrons do not annihilate on the container walls 
and are easy to produce, the final neutral plasma temperature is tuned by 
setting the electron beam temperature. 

We assume that the positrons are produced continuously through the inter- 
action of and intense beam of 10 MeV electrons with a thin tungsten foil [20] 
with a rate, after solid Neon moderation, of 10^^ s~^. The positrons are cooled 
to room temperature (~ 25 meV) in a buffer gas. In the slow loading mode 
the positrons are taken at the exit of the buffer gas section while in the fast 
loading mode they are first accumulated in a Greaves-Surko trap [16,14,15] 
where they are stored and cooled to ~ 2 meV. The low temperature of the 
positrons is a feature of the Greaves-Surko trap: the strong magnetic field 
is produced by a supraconducting magnet and the positrons are in thermal 
equilibrium with the magnet container at a temperature of 2 meV. Such a low 
temperature enhances the 3 body reaction which absorbs the positrons if the 
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Fig. 2. The MCEO trap electric potential (a) and magnetic flux (b) set for antihy- 
drogen synthesis [17]. 

neutral plasma were to be stored in this trap. 

The preferred device for accumulating the e"^ in a neutral plasma is the MCEO 
trap [17] which does not require a very low temperature such as the one 
reached in a Greaves-Surko trap. Therefore the heating of the positrons can be 
accomodated and is beneficial to accumulate enough positrons in the neutral 
plasma. 

The positrons are guided by a 100 Gauss magnetic field, from the Greaves- 
Surko trap (or the buffer gas), to the neutral plasma. The injection pipe is 
shielded to minimize its leaking field. 

The original MCEO trap is a magnetic cusp with an octupolar electric field. 
It is made of several pairs of coils of 10 cm radius arranged symmetrically 
with respect to the center of the trap. The most central coils are located at 4 
cm from the center. The other coils are separated from each other by 5 cm. 
The currents in each pair of symmetric coils are of opposite sign. The current 
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amplitude can be varied for each pair. 

It is designed to mix both positrons and antiprotons. It is to be implemented 
in the ASACUSA beam line at CERN under the name of Musashi trap. In 
our application this trap is meant to mix positrons with electrons. Therefore 
a scaled down version has to be designed. 

The scaled down version of the MCEO trap has the following parameters: the 
magnetic coils radius is 1.5 cm, the most central coils are at 1.2 cm from the 
center, the next coils are separated by 1.5 cm, the total current for each set of 
coils is 4 k Ampere-turns, the octupolar electrostatic field is set by a central 
electrode at 70 V and a length of 1.0 cm. The key feature of this trap is to have 
a null field region in its center (figure 2). Once a charged plasma is loaded, 
it acts as a potential well for the opposite charge. We use electrons to create 
the well for the positrons. The neutral plasma will be located at the center of 
the well. The e~ and beams are injected along the x axis. The octupolar 
electrostatic field of the MCEO trap is simulated using an analytical form [21]. 

A simple trap made of two coils acting as mirrors [22] is also considered but 
the injection of the positrons is more difficult. 

In order to extract the positrons from the neutral plasma and accelerate them 
toward the crystal, an electrostatic field is used. The electrode configuration 
is shown in figure 3. 

The electrodes are arranged in two parallel planes. The cathode is made of 
Al cristal. Its surface is 100 /im x 1 cm. The neutral plasma is located in the 
volume between the two planes. This configuration will be referred as Charge 
Focusing Aluminium Converter or CFAC. 

This configuration is inspired by a technology developped for high energy 
physics detectors called the Micro Strip Gas Chamber [23], or MSGC. The 
distances between electrodes and field magnitude are taken from this devel- 
opment. 

In the CFAC, the electrostatic field of the MSGC is modified by introducing 
4 wires to smoothe the side field as shown in figure 4. 

The CFAC was put in the center of the trap by the simulation described below 
in order to get an estimation of the positronium density and time distribu- 
tion. The insertion of metallic conductors in the center of the MCEO trap is 
unrealistic. An integrated design is proposed in the discussion section. 

A fortran program was written to simulate: 

• the stability of the neutral plasma during the 10 /iS loading time. 
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Fig. 3. The CFAC electrodes geometry. 
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Fig. 4. The CFAC electrostatic potential: (0,0) is the center of the crystal surface. 



• the acceleration by the CFAC electrostatic field and the time distribution 
of the positrons hitting the Al crystal. 
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Fig. 5. Positron tracks from the plasma to the CFAC in cm. 
2.1 The target simulation 

In the MCEO trap the positrons are injected along its axis, parallel to the 
magnetic field which guides them. 



2.1.1 The CPA C electrostatic field 

The CFAC electrostatic field is computed using analytical formulae [23]. The 
gap of the CFAC, i.e. distance between the drift electrode and the cathode 
planes is 1.5 cm. The cathode plane is located at 1 cm from the x axis and 
parallel to it. The drift electrode size is 1 cm x 7 mm. The cathode, made of 
the Al(lll) crystal has a 100 /xm width. The distance between the edge of 
the cathode and the next electrode in the cathode plane is 200 /xm. The drift 
electrode plane being at 1600 V, the cathode is grounded and the bias voltage 
is 1550 V. The plasma temperature is 30 meV. 

During the loading phase the CFAC field is set "off" . The acceleration phase 
starts when the CFAC electrostatic field is set "on" while the octuplole elec- 
trostatic field and the cusp magnetic field are set " off" . 

The positron capture and the emission of positronium atoms is a fast process, 
the delay of which can be neglected at the nanosecond scale [24]. Therefore 
the time distribution of the positrons hitting the crystal is also the time dis- 
tribution of the positronium atoms exiting the crystal surface (figures 5,6). 
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Fig. 6. Arrival time distribution of the positrons on the crystal in nanoseconds. 

The positronium atoms exiting the crystal are in the ground state. The mag- 
netic field couples the singlet and triplet states [25], reducing the Positronium 
life time to a minimum of 0.125 ns. 

In an experimental test with a MCEO trap and with electrons the radial field 
on the central electrode surface was 360 Gauss [26]. This trap is designed to 
hold both positrons and antiprotons. When antiprotons are used the magnetic 
field is much stronger but in the beam to cloud configuration the antiprotons 
cross the trap axis and are not kept. Hence with a weak field (~ 3.5 Gauss) 
we take as effective lifetime 0.200 ns. This life time corresponds to a travel 
distance of 200 fim. The target volume is defined as a box parallel to the 
crystal, which starts at a distance of 100 /im. The target is 100 /im large, 1 cm 
long and 100 /im thick, which gives a volume of 10~^ cm^. This target volume 
is centered at a distance of 150 /im from the crystal surface. 

An estimate of the target density is obtained by assuming all Ps atoms to fiy 
perpendicularly from the crystal surface [11]. The number of positrons stored 
in the plasma is A^e+ = 10^^. The positron fiux from the plasma to the crystal 
has a duration of ^ 10 ns. The Ps atoms spend only 0.100 ns in the target 
volume, hence the target density is given by: 



The target volume could also be counted starting from the crystal surface with 
some antiproton losses (antiprotons hitting the crystal) . 



10^1 0.100 
10^ 10 



^Ps ^CFAC = 0.25 10 



cm 



-3 



(1) 
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The simulation of the CFAC gives an efficiency of ecFAC = 0.5 for a positron to 
reach the crystal. The positronium emission efficiency from the crystal surface 
is taken as ep^ = 0.5. 

A correcting factor due to the positronium decay between the crystal and the 
target volume will reduce the density by less than a factor 2. This correction 
depends on the magnetic field strength. 



2.2 Positron losses 



When mixing e~ and e+ with a density above ^10^ cm~^ and a tempera- 
ture of a few meV, the plasma collapses into positronium through a 3 body 
interaction. 



e- + e+ + e± ^ Ps(nT)+e^ (2) 



The rate of other reactions, namely the radiative recombination and the di- 
rect annihilation are negligible. This reaction and the subsequent positronium 
annihilation destroy the positrons. 

The number of positrons in the neutral plasma as a function of time is given 
by 1; 



TVe = X tanh(\/aA t) (3) 
V A 

limiV.W = ^ (4) 



Where a is the injection rate and A is a constant related to the 3 body reaction. 

The tables obtained in the appendix show that in the worst case, in ^ 2 s and 
at a temperature equal or above 4 meV it is possible to accumulate 2 10^^ in 
the plasma. For less than 10 seconds of accumulation there is little variation 
with temperature above 4 meV. Once the injection of positrons ends and 
during ^ Is the loss of positrons is small and on the time scale of a few /iS it 
can be neglected. 

In order to have a negligible loss by annihilation on remaining gas a vacuum 
environment is required with a residual gas partial pressure below ^ 10~^ torr. 



^ see computation in appendix 
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Fig. 7. Hydrogen formation cross section in Ps + p ^ H + e+ collisions [27] . 
Closed circles represent the calculation of Yamanaka and Kino; lines, other 
models referenced in the same paper; crosses, the experiment of Morrison et 
al. [28,31]. The abscissa is the energy of the proton in the center of mass frame. 

3 Antihydrogen production 

The charge exchange on positronium to produce the atom is a resonant re- 
action. The cross-section calculations [27,29,30] have been confirmed by an 
experiment [28,31]. 

Once the H atom is produced, it captures a positron to make the ion. There 
is a non resonant channel, the two body recombination, and a resonant one, 
the charge exchange on positronium. 

The choice of the antiproton kinetic energy is given by the overlap of the 
resonant cross sections. In the region of interest, the antiprotons have a kinetic 
energy of ^ 13 keV (figure 7). The atom production cross section obtained 
experimentally [28] is 7.8 Tva^ ^ . At this kinetic energy, the computed cross 
section [32] for the ion production is ^ 0.05 iral (figure 8). 

The (anti) hydrogen production rate normalized to the (anti)proton flux is: 

^Ps X cr^^^P (5) 

The (anti)proton speed at 13 keV is vp — 1.58 lO^cms"^. Therefore the nor- 
malized rate is 0.28 10^ s~^. The target length is 1 cm, hence the crossing time 

^ Here ao is the Bohr radius and ttOq = 0.880 x 10"-^^ cm^. 
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Fig. 8. H~ production cross section in H(n = 1) + Ps(^' = 1) — > H~ + e"^ colli- 
sions [32]. The lines are two variants of a two-channel calculation by Biswas. The 
abscissa is the energy of the Ps in the center of mass frame. 

by the (anti) proton beam is 6.3 ns and the number of antihydrogen atoms pro- 
duced by each antiproton crossing the target is 0.0017. 

The (anti) hydrogen ion production rate normalized to the antihydrogen flux 
is: 



(6) 



The normalized rate is 0.18 10^ s~^. After convolution for the H production 
in the same 1 cm long target, the number of ions produced for each antiproton 
crossing the target is 10~^. 

At the edge of the crystal and at a distance of 150 fim the CFAC electro- 
static field component transverse to the antiproton path is ^ 0.12 Mega 
Volts/m. Therefore an antiproton aiming at the edge of the crystal will be 
pulled transversally by 240 /xm on a distance of 1 cm. The transversal electro- 
static field which is null on the axis and grows toward the edges of the target, 
will cause a maximum pull ^ 0.5 mm and a beam divergence of ^ 50 mrd. 
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4 Discussion 



In the integration of the devices (CFAC, MCEO trap, buffer-gas, Greaves- 
Surko trap) into a single experiment there is room for changing their pa- 
rameters: gains from tuning the device parameters and losses are difficult to 
simulate and have to be established experimentally. When changing the pa- 
rameters in the simulation these gains/losses are within a factor 10. 

Still there are many technical questions which are not addressed by the sim- 
ulation. We shall now hst some of them. 

The field in the buffer gas trap is 1.5 kGauss and 5 T in the Greaves-Surko 
trap: a set of iron plates and auxilliary fields shall be designed to reach the 
100 Gauss guiding solenoidal field [16]. 

Positron plasmas with a density of 10^°e+cm~^ at meV temperatures have been 
stored in a Greaves-Surko trap with a few 10^ trapped e^. The next generation 
of multicell traps [16,14,15] is expected to store and cool to meV temperature 
lO^^e"^. The fast loading mode cannot be implemented without traps that can 
store at least 10^^ In this fast mode, once trapped, the positrons shall be 
extracted in a single pulse. This step has to be experimentally established. 

Such traps cannot be filled with e+ emitted by a radioactive source in less 
than a few weeks. An intense source of slow positrons such as the one we 
proposed [20], using an intense 10 MeV electron beam impinging at a small 
incidence angle on a thin tungsten foil is expected to fill the trap at a rate of 
lO^^e+s-^ 

In the slow loading mode, no trap is needed and only a buffer gas is required. 

At high density, the neutral plasma is not transparent to the electro-magnetic 
field. The extraction of the positrons by the CFAC field in ^ 10 ns requires a 
pulse shaping technique [33] and shall be tested experimentally. 

In order to implement a gravitation experiment with a scheme such as the 
one discussed by J. Walz and T. W. Hansch [7] constraints have to be met by 
the antiproton source. The ions are produced in a beam at 13 KeV with an 
energy and angular spread dominated by the initial (anti) proton energy spread 
and the ^ 50 mrd divergence due to the CFAC field. The angular spread due 
to the charge exchange reactions is negligible. The angular spread shall be 
taken into account to decelerate and collect the into an ion trap with high 
efficiency [34]. Antiproton traps already exist, but a trap with lO^p capable 
of producing a beam of 13 keV ± 1 keV energy with a divergence lower than 
0.01 rd (ie 100 /im on 1 cm path) is still to be demonstrated. 



13 



Fig. 9. Modified central electrode of a MCEO trap to produce an approximate CFAC 
field. 

For antihydrogen creation, the very large plasma densities discussed above 
allow to transform the antiprotons with a relatively high efficiency. But an- 
tiprotons are difficult to produce. However, in order to test the process one can 
use protons and an existing trap: the lower plasma density being compensated 
by a larger number of protons (say 10^ p per beam pulse). 

Another aspect of this path for a gravity experiment is that the reactions 
involved are at all steps charge symmetric: by switching from protons to an- 
tiprotons on the time scale of 1 hour, it is possible to measure g in the same 
gravitational field with little or no tide effect due the movement of the moon 
and other masses. 

The integration of the CFAC into a MCEO trap requires some modification. 
We foresee to cut the central electrode in 4 (or more) sections separated by a 
thin insulator (figure 9). In the trapping mode all the sections are set at the 
same potential acting as a single electrode. In the CFAC mode the potentials 
are set to produce an electrical field similar to the one in figure 3. 

The ion beam is needed for the gravity experiment, but a large number 
of H atoms are also produced in the beam. Using a laser to stimulate the 
transition toward the n = 2 level one gets a beam in the ^10-20 keV range. 
Then the Separated Oscillatory Field method [35,36,37] allows to measure the 
25^1/2 — 2P3/2 fine structure. 

Appendix 



14 



The technique described in this article can be used for experiments on positro- 
nium BEC and a 511 keV laser because it could reach a positronium density 
of ^ 10^^ or ^ 10^^ cm~^. The appendix gives a numerical estimation of the 
density as a function of time and temperature. 

The maximum number of positrons available to interact with an Aluminium 
crystal is the result of the equilibrium between the injection rate and the 3 
body reaction rate. This reaction (see equation 2) and the subsequent positro- 
nium annihilation destroy the positrons. 

The time evolution of the number of positrons in the plasma is given by the 
injection rate and the 3 body reaction rate [38]: dN^ = a dt — dt 

where N^. is the number of in the plasma volume, a is the injection rate 
and rs is the 3 body reaction rate given by: — A x nl x Ve x 



where the electron thermal velocity = \/ and the impact parameter 



b = ff are numerically: Ve{cm.s ^) = 4.19 10^ JT{eV) and 6(cm.eV ^) = 



1.44 10-7T(eV), resulting in rsis'^) = 2.6 10"^^ A x nl / T^'^ 

where A is a parameter which varies with the magnetic field, rie is the density 
in cm~^ and T is the plasma temperature in eV. 

The neutral plasma volume is ^ 30 mm^. 

Let V be the plasma volume and let's define A as A = 2.6 10"^^ A/{V'^ T^-^), 
so that r3 = A N^. 

With the initial condition Ne{t = 0) = 0, we get: dNe = (a - A A^^^) dt 

Which leads to: N^^ ^x tanh(\/aT t) , lim^^oo N^it) = 

For the two extreme situations, B — Q and B — oo^ the A parameter is 
respectively equal to 0.76 and 0.07. 

In the computations the ion was supposed to be much heavier than the elec- 
tron and therefore its trajectory in the field was neglected [38]. Recently a 
computation was made with a proton taking into account its trajectory [39] 
for fields of a few Tesla: the variation of the parameter A compared to the 
infinite field value was less than a factor 2. Here the ion is the positron which 
is much lighter and the magnetic field is a few Gauss only. Therefore the 
positron behaves as a heavier ion in a stronger field. The expected behavior 
is then an intermediate one, between the above computation with field and 
with infinite field. 

A computation involving only electrons and positrons predicts a three body 
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Table 1 



Number of 6+ accumulated in the 30 mm^ neutral plasma with an injection rate 
a = 10"s-^ for A = 0.76 (B = 0) 



T (meV) 


2 


4 


10 


25 


50 


100 


0.1 s 


9.99 10^ 


9.99 10^ 


9.99 10^ 


9.99 10^ 


1.0 10^° 


1.0 10^° 


1 s 


9.91 10^° 


1.0 10" 


1.0 10" 


1.0 10" 


1.0 10" 


1.0 10" 


2 s 


1.93 10" 


2.0 10" 


2.0 10" 


2.0 10" 


2.0 10" 


2.0 10" 



Table 2 

Number of e+, accumulated in the 30 mm^ neutral plasma with an injection rate 
a = 10^^5-1 for A = 0.07 (B = oc) 



T (meV) 


2 


4 


10 


25 


50 


100 


A = 0.76 
A = 0.07 


1.81 10" 
5.95 10" 


8.60 10" 
2.83 10^2 


6.76 10^2 
2.23 10^^ 


5.31 10^^ 
1.75 10^^ 


2.53 10^^ 
8.32 10^^ 


1.20 10^^ 
3.96 10^5 



Table 3 



= maximum number of e^, that can be accumulated in the 30 mm^ neutral 
plasma with an injection rate a = 10^^ for A = 0.76 (B = 0) and A = 0.07 (B 
= oc) 

recombination rate twice greater than for protons and electrons [40] . 



The computed number of positrons accumulated inside the neutral plasma for 
several temperatures is tabulated for both values of A. 

When the injection in the plasma ends, the evolution is given by the rate of 
the 3 body reaction: 

dNe = -rs dt = -A Nl dt 
M = ^io 

i+7Vo A t 
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